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Abstract. 

The Bohr-Sommerfeld approximation to the eigenvalues of a one-dimensional 
quantum Hamiltonian is derived through order (i-e., including the first 
correction term beyond the usual result) by means of the Moyal star product. 
The Hamiltonian need only have a Weyl transform (or symbol) that is a 
power series in h, starting with fp , with a generic fixed point in phase space. 
The Hamiltonian is not restricted to the kinetic-plus-potential form. The 
method involves transforming the Hamiltonian to a normal form, in which it 
becomes a function of the harmonic oscillator Hamiltonian. Diagrammatic and 
other techniques with potential applications to other normal form problems are 
presented for manipulating higher order terms in the Moyal series. 



PACS numbers: 03.65.Sq, 02.40.Gh, 03.65.Vf, 02.20.Sv 
1. Introduction 

In this article we use the Moyal star product to derive the Bohr-Sommerfeld 
approximation for the eigenvalues of the bound states of a one-dimensional 
Hamiltonian, including higher order corrections in h. We begin by quoting the result. 



En — 



,(1) 

= (n+l/2)h 

which uses the following notation. £'„ is the n-th eigenvalue of the quantum 
Hamiltonian H, which has Weyl transform (or "symbol") H. The latter is treated 
as a classical Hamiltonian with action-angle variables (A, 0), and is regarded as a 
function of the action A. The frequency of the classical motion is lo{A) = dH/dA, 
and the notation {H, H}2 refers to the second Moyal bracket, defined in (|A.4c|l below. 
This Moyal bracket is otherwise twice the Hessian determinant of the Hamiltonian, 

{H, H}2 = 2[H^xxH^pp — (H^xp)^)]- (2) 

The angle brackets {■ ■ ■)if> represent an average over the angle 0. This result is discussed 
further in subsection 15. II 

This paper assumes some background in the Wigner-Weyl formalism and the 
Moyal star product. A sampling of references in this area includes Weyl (1927), 
Wigner (1932), Groenewold (1946), Moyal (1949), Berry (1977), Balazs and Jennings 
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(1984), Hillery et al (1984), Littlejohn (1986), McDonald (1988), Estrada et al (1989) 
and Gracia-Bondia and Varilly (1995). 

Our strategy is to use the Wigner-Weyl symbol correspondence and the series 
representation of the Moyal star product to transform a given Hamiltonian into a 
function of the harmonic oscillator, at least in a neighborhood in phase space of a 
fixed point of the classical Hamiltonian. Notable aspects of this calculation are the 
higher order terms in the Bohr-Sommerfeld formula, the use of the Moyal star product 
to achieve a normal form transformation, and the development of diagrammatic and 
other techniques for handling higher order terms in the Moyal series. 

In one-dimensional problems with Hamiltonians of the kinetic-plus-potential form, 
it is straightforward to extend standard WKB methods to higher order in h to find 
corrections to the usual Bohr-Sommerfeld formula. Results may be found in Bender 
and Orszag (1978). Several other approaches to the calculation of higher order terms 
have been proposed in the literature (Maslov and Fcdoriuk, 1981; Voros, 1977, 1989; 
Kurchan et al 1989). Our approach is characterized by the use of the Weyl symbol 
correspondence for the representation of operators, and the Moyal star product for the 
manipulation of those operators. In a subsequent paper we shall extend our methods 
to the case of multidimensional, integrable systems (torus quantization). 

Recently Colin de Verdiere (2004) has presented another approach for calculating 
higher order corrections to the one-dimensional Bohr-Sommcrfeld formula that is 
based on the Weyl symbol correspondence and the Moyal star product. Colin de 
Verdiere concentrates on the spectrum of the operator, which is manipulated by 
representing traces of operators as integrals over phase space. In his approach there is 
no need to examine eigenfunctions. From an algorithmic or computational standpoint, 
his method is quite simple, the simplest we have seen for calculating higher order 
terms (certainly simpler than ours). On the other hand, our approach does provide 
explicit representations of the transformations necessary to construct eigenfunctions. 
Another difference is that our method can be generalized to the multidimensional case, 
whereas we could not see how to do that with the method of Colin de Verdiere. The 
generalization of our method was not trivial, however, and it may be that with more 
effort his could be generalized, too. 

The use of the Weyl symbol correspondence for representing operators means 
that one can handle a wider class of Hamiltonians than the kinetic-plus-potential 
type (second order differential operators). For example, integral as well as differential 
operators are allowed. This is important in applications, such as in plasma physics, 
where such operators arise, and also for the multidimensional generalization where not 
all the operators of the commuting set need be second order differential operators. 

More importantly, the use of the Weyl symbol correspondence means that 
the calculations take place in phase space, not configuration space. For example, 
there is no trouble with caustics or turning points. (Phase space methods are 
not necessarily free of caustic difficulties, but caustics make no appearance in 
our approach.) For another example, the usual (lowest order) Bohr-Sommerfeld 
formula has an obvious invariance under arbitrary canonical transformations, since 
the energy eigenvalues are expressed in terms of the area of the classical orbits 
in phase space. To make this statement precise, however, one would have to 
say precisely what the classical Hamiltonian corresponding to a given quantum 
Hamiltonian is, over a wider class of operators than those of the kinetic-phis-potential 
type (a form that is not invariant under canonical transformations). The Weyl 
symbol correspondence does this, and provides a formalism that is covariant under 
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linear canonical transformations. Thus, in approaches based on the Weyl symbol 
correspondence, the lowest order Bohr-Sommerfeld energy eigenvalues are invariant 
under aribitrary canonical transformations, and the next and all higher order terms 
in the H series are invariant under linear canonical transformations. 

The basic idea of this paper arises from the usual, lowest order Bohr-Sommerfeld 
formula, which states that the eigenvalues of a quantum Hamiltonian are given 
approximately by setting A = {n + 1/2) ft in the classical formula expressing the 
classical Hamiltonian _ff as a function of its action A, H = f{A). This formula 
suggests that the quantum Hamiltonian is a function of a quantum "action operator," 
of which the classical formula is a lowest order representation by means of symbols, 
and that the eigenvalues of the action operator are (n+ l/2)fi. Since these are also the 
eigenvalues of the harmonic oscillator (of unit frequency), the suggestion is that the 
action operator is unitarily equivalent to the harmonic oscillator Hamiltonian. If this is 
so, then the quantum Hamiltonian is unitarily equivalent to a function of the harmonic 
oscillator. In this paper we find that these suggestions are borne out, and we are able 
to construct explicitly the unitary operator (or at least the symbol of its generator), in 
the sense of a formal series in h, which transforms a given quantum Hamiltonian into 
a function of the harmonic oscillator. We only require that the quantum Hamiltonian 
have a "slowly varying" (defined below) Weyl symbol, and that the symbol have a 
generic extremum (fixed point) at some point in phase space. The classical analog of 
the unitary transformation we construct is a canonical transformation that maps the 
level sets of the classical Hamiltonian around the extremum (which are topological 
circles) into exact circles about the origin. The latter, of course, are the level sets of 
the harmonic oscillator. 

The transformed Hamiltonian can be regarded as a normal form, that is, a 
standard or especially simple form achieved by means of a transformation. In this 
case, the normal form is a function of the harmonic oscillator, achieved on the level 
of 7i-series for the Weyl symbol and brought about by unitary transformations. The 
class of Hamiltonians which can be brought into this normal form are those whose 
symbol has certain properties, notably that of having a family of level sets of circular 
topology surrounding a generic fixed point. This is the most generic case for bound 
states in a one-dimensional system, and therefore the one to examine first. But one 
can imagine that there are other normal forms that apply in other cases, for example, 
if the fixed point is not generic perhaps a standard quartic oscillator is a normal form. 
Or if in a certain region of phase space the symbol has a separatrix surrounding two 
islands of equal area and crossing in one unstable fixed point, perhaps a standard, 
symmetric double well oscillator is a normal form. Certainly at lowest order in h the 
unitary transformation is represented by a canonical transformation that preserves 
area in phase space, so the separatrix of the Hamiltonian and that of the normal form 
must enclose the same area if they are to be unitarily equivalent. Whether this is 
enough to guarantee the equivalence of the operators (that is, the equality of their 
symbols to all orders in h) is an open question, as far as we know. 

Normal form transformations of operators or sets of operators, either by star 
product methods or by Fourier integral operators, have been around for some time, in 
fact we ourselves have used the star product to diagonalize or block diagonalize systems 
of coupled wave equations (Littlejohn and Flynn, 1991) and to study problems of 
mode conversion (Littlejohn and Flynn, 1992, 1993). Such methods are well suited to 
the development of systematic perturbation methods for quantum adiabatic systems, 
such as molecules in the Born-Oppenheimer approximation (Littlejohn and Weigert, 
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1993). Similar normal form transformations for coupled wave equations have also been 
discussed by Braam and Duistcrmaat (1995), although with the idea of using Fourier 
integral operators to carry out the transformations instead of star products. The star 
product approach to normal form transformations for coupled wave equations (WKB 
on vector bundles) was put on firmer mathematical foundations and generalized in 
various ways by Emmrich and Weinstein (1996) and Emmrich and Romer (1998). 
More recently, Colin de Verdiere et al (1999) have studied normal forms for mode 
conversion (generalized Landau-Zencr transitions) and Colin de Verdiere and Parisse 
(1999) have studied them for problems involving separatrices. We also expect normal 
form theory to have applications in transition state theory (Creagh, 2004). 

In recent years there has arisen the subject of deformation quantization, 
which involves generalizations of the Moyal star product to nontrivial phase spaces 
(symplectic or Poisson manifolds). The phase space used in this paper, upon 
which the Wcyl symbol correspondence is based, is considered trivial. The general 
idea is to deform the commutative algebra of multiplication of functions on the phase 
space into a noncommutative but associative algebra, where h is the deformation 
parameter and where the new multiplication rule is the generalized star product. It is 
also required that the order h term in the symbol representation of the commutator 
be proportional to the Poisson bracket. The new algebra is then interpreted as an 
algebra of operators on a quantum system, the "quantized" version of the classical 
phase space. In one approach, the star product is represented as a formal power series 
in h, a generalization of the Moyal formula, and one must work out the terms of 
the series subject to the constraint of associativity and the appearance of the Poisson 
bracket at first order. Basic references in this area include Bayen et al (1978), Fedosov 
(1994) and Kontsevich (2003). These works show an interesting geometrical structure 
associated with the higher order terms in the Moyal star product, which has stimulated 
our interest in higher order terms in the Bohr-Sommerfeld formula. Our results may 
be relevant for attempts to understand eigenfunctions of an operator geometrically as 
noncommutative versions of tori. 

It should be possible to read the main body of this paper, skipping the appendices, 
to obtain an overview of our calculation. The appendices, however, are needed for the 
details, including notational conventions. 

2. The setup 

Let if be a Hermitian operator (the "Hamiltonian" ) in a one- dimensional quantum 
system, that is, H acts on wave functions tl}{x), a; e M (the Hilbert space is L^(M)). 

We uniformly use hats ( ) over a letter to denote operators, whereas a letter without 
a hat represents the Weyl transform (or Weyl "symbol" ) of the operator. For example, 

H{x,p) = J ds e-'P'^'' {x + s/2\H\x - s/2) (3) 

and 

H = j dj^^,i.i.'-.)inH \x){x'\ (4) 

illustrate the Weyl transform and its inverse in the case of the Hamiltonian. We regard 

H as the "classical Hamiltonian," defined on the phase space K.^, with coordinates 
{x,p). We denote these coordinates collectively by z'^ = {x,p), = 1,2. 
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We assume that the symbol H has an expansion in h beginning with the power 

H = Ho + hHi+ tfH2 + ..., (5) 

where each Hn is independent of h. An operator whose symbol possesses this type 
of expansion will be called "slowly varying". Not all operators are slowly varying; 
for example, the unitary operator exp{—iHt/h) is not. The leading term {Hq in the 
example above) of the symbol of a slowly varying operator will be called the "principal 
symbol." 

We assume H is smooth and has a generic extremum (a fixed point) at some 
point of phase space. The fixed point need not be at p = 0, nor docs H need to be 
invariant under time-reversal {p —p). An extremum is considered generic if the 
Hessian matrix H^^i, of the Hamiltonian is nonsingular at the extremum. Here and 
below we use comma notation for derivatives, for example. 

For example, the fixed point (x, p) = (0, 0) of the quartic oscillator {V{x) = x^) is not 
generic, because the Hessian matrix has rank 1 at the fixed point. 

It is convenient in what follows to assume that the extremum is a minimum 
(always the case for kinetic-plus-potential Hamiltonians). If not, we replace H by —H 
at the beginning of the calculation. 

Radial equations (on which x is the radial variable r > 0) are excluded from 
our formalism, because the Weyl symbol correspondence is not defined in the usual 
way on the half line, and because the centrifugal potential is singular. We believe 
the best way to handle such problems within a formalism like the one presented in 
this paper is by reduction from a problem on a higher dimensional configuration 
space R" under some symmetry, typically SO{n). Nor are singular potentials such 
as the Coulomb potential covered by this formalism, because singularities generally 
invalidate the Moyal star product expansion in h, itself an asymptotic expansion. The 
usual lowest order Bohr-Sommerfeld formula usually does give correct answers for 
singular potentials, at least to leading order in h, but the structure of the higher order 
terms (which powers of H occur, whether the corrections can be represented by powers 
of h at all, etc.) presumably depend on the nature of the singularity. 

In view of our assumptions, the classical Hamiltonian H has level sets in some 
neighborhood of the fixed point that are topological circles. We concentrate on this 
region of phase space, and ignore any scparatrices and changes in the topology of the 
level sets of H which may be encountered further away from the fixed point. 

For convenience we perform a canonical scaling on the coordinates {x,p) (or 
operators {x,p)) to cause them both to have units of action^^^. For example, in the 
case of the ordinary harmonic oscillator, we would write x' = y/nwjx, p' = p/y/muj, 
and then drop the primes. 

The strategy will be to perform a sequence of imitary operations that transform 
the original Hamiltonian H into a new Hamiltonian that is a function of the harmonic 
oscillator Hamiltonian, at least in the "microlocal" sense of the symbols in the 
neighborhood of the fixed point. Since unitary transformations do not change the 
spectrum of an operator, the new Hamiltonian will have the same spectrum as the 
old one. But since the new Hamiltonian is a function of the harmonic oscillator 
Hamiltonian, its eigenvalues are easy to compute. 
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The transformations will proceed in two stages. In the first stage, we perform 
a "preparatory" transformation that maps H into a new Hamiltonian K that is a 
function of the harmonic oscillator Hamiltonian at lowest order in h. We follow this 
by a sequence of near-identity unitary transformations that transform K into a new 
Hamiltonian M that is a function of the harmonic oscillator Hamiltonian to all higher 
orders in h, at least formally. Thus, the stages are 

H k ^ M. (7) 

What we mean by the harmonic oscillator Hamiltonian is really the action of the 
harmonic oscillator, given in operator and symbol form by 

i=l{i'+f), (8) 

It turns out that an operator is a function / if and only if its symbol is a function 
of J, as will be discussed more fully below, although the two functions are not the 
same beyond lowest order in H. Thus, to ensure that the transformed Hamiltonian is 
a function of /, we require that its symbol be a function of /. 



3. The preparatory transformation 

The preparatory transformation (the first arrow in |(7J)) is the most difficult, because 
it is not a near-identity transformation and cannot be handled by Lie algebraic (power 
series) methods. This transformation will transform H into another Hamiltonian k 
whose symbol is a function of / plus terms of order and higher. Thus, the principal 
symbol of k will be a function of /. The preparatory transformation only makes the 
leading order term in the symbol of a function of /, not the higher order terms. 

3.1. Imbedding H and U in a family 

Let H be given. Ultimately, we shall seek a unitary transformation U such that the 
new Hamiltonian K, defined by 

k = UHU\ H = u''ku, (9) 

has a symbol K that is a function of / plus terms of order ft^ and higher. 

For the moment, however, it is conceptually simpler to imagine that H and U are 
given, and to seek a means based on Weyl symbols of computing k, without regard 
to the functional form of K. We do this by imbedding [/ in a family, < e < 1, that 
is by assuming that there exists a smooth family of unitary operators ?7e, such that 

The family C/e can be seen as a path in the group of unitary transformations that 
act on Hilbert space, connecting the identity and the final U. We do not assume e is 
small, and we do not carry out any power series expansion in e. We imbed H in a, 
similar family, defining 

H.^iilku,, (11) 

so that 

i7, = ( ^ ^^' = 0' (12) 

1 iJ if e = 1. ^ ' 
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Table 1. Notation for operators, symbols and functions depending on e. 



e = 1 K Kn Id'' e I 
any e Ue H„c Af^ 

e = l U H 5„oH Z^' (p A 



One might say that the e-evolution runs backwards, since K evolves into H as e goes 
from to 1. As always, Weyl symbols of the operators above are denoted without the 
hat, for example, C/, if, U^^ H^. There are several operators, symbols and functions 
in this paper that depend on e, the notation for which is summarized in Tabled We 
shall be interested in calculating H^^ from which K follows by setting e = 0. 

We obtain a differential equation for by differentiating which gives 

^ = ^[Ge,ii.], (13) 

where the Hermitian operator Ge (the "generator" ) is defined by 

G, ^ihUl^^~ih^U,^G\. (14) 

We assume that is slowly varying. We shall solve (|13|l by converting operators to 
symbols and using the Moyal product formula. See [Appendix A| for the Moyal star 
product and the diagrammatic notation we shall use for the functions and operations 
that arise from it. 

3. 2. Differential equations for H^^ and Hne 

We now transcribe H13|l to symbols and substitute (|A.6p . This gives a differential 
equation for the symbol H^, 

^ = -{G„ i/J + g{Ge, H^h ~ i?j5 + . . . , (15) 

which is subject to the boundary condition = H a.t e = 1. We express the solution 
of this equation in terms of a set of new functions or symbols, Hq^, i?2c, etc., which 
are required to satisfy the differential equations, 

^-{i/o.,GJ = 0, (16a) 
de 

^-{H2,,G,} = ^{G,,Hoe}3, (16&) 

de 24 
^ - {ff4.,GJ = l{G,,H2.h - jko^G^Ho^h, (16c) 
etc., and the boundary condition Hm = Sno H at e = 1. Then we have 

H, = Hoe + h^H2e + h-^H^e + ■■■■ (17) 

This is not an expansion of in powers of fi. as in jSJ, because the functions Hne are 
themselves allowed to have a dependence on h. But each of these is slowly varying, 
so that if the series (I17|l is truncated, the order of the omitted term is given by the H 
coefficient. Finally, we define Kn = evaluated at e = (see Table so that we 
have an expansion of the symbol K of K, 

K = Kq + h^K2 + h'^K^ + . . . . (18) 
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The solutions of (|16flfl - Hlt)c|l can be expressed in terms of a certain e-dependent, 
classical canonical transformation, z''^{e) — Z^{z)^ where z and z' are the old and 
new variables, and is the transformation function. The family of canonical 
transformations reduces to the identity at e = 0, while at e = 1 we shall denote the 
transformation simply by (without the e). See Table The transformation 
will be defined momentarily, but it turns out that the left hand sides of (|16af) - (|16cll 
are convective derivatives along the associated Hamiltonian flow. Equation l|lt)a|l is 
a homogeneous equation for the unknown i?oej and the others are inhomogeneous 
equations with driving terms determined by lower order solutions. The structure of 
the system is that of a Dyson expansion, in which the canonical transformation Z^ 
specifies a kind of interaction representation. The definition of Z^' requires some 
notational understandings that are presented in [Appendix B| 



3.3. The canonical transformations Z and Z^ 

The canonical transformation is defined as the solution of the functional differential 
equation, 

dZ^^ 

^ = {^r,Gj, (19) 

subject to the initial conditions, Z^ = Id^ at e = 0, and we define Z'' = Zf at e = 1 
(see Table The functions Zf so defined constitute a canonical transformation, for 
if we compute the e-derivative of their Poisson brackets among themselves, we find 

^{Zf, Zn = {{Zif, GJ, zn+{Zif, {Z^, GJ} = -{G„ {Zf , Zn}, (20) 

where we have used the Jacobi identity. These are subject to the initial conditions 
{Zf , Z^} = J^'' at e = 0. But since J^" — const, the initial conditions are the solution 
for all e, as shown by direct substitution. 

The canonical transformation Zf is not generated by G^ regarded as an e- 
dependent Hamiltonian function, but rather by G^ = Ge o Z~^. That is, if we write 
z^(e) = Z^(zo) for the solution of Hamilton's equations, 

^ = J^''G;,(z), (21) 

then the functions Zf satisfy 
dZf^ 

^ = J^'' G;, o Z, = {Id^ G^} o Z, - {Z,^ GJ, (22) 

which agrees with (|19l) . In the final step we have used an important property of the 
Poisson bracket, namely, that if A and B are any two functions and Z is a canonical 
transformation (symplectic map), then 

{A,B}o Z ^ {Ao Z,Bo Z}. (23) 



3.4. Notation for e-derivatives 



The following notation will be useful for carrying out differentiations and integrations 
in the interaction representation, specified by composing a function with Z~^. 
For any function on phase space, possibly e-dependent, we define 



De 



de 



(FeOZr^) 



O Zf 



(24) 
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for a kind of derivative operator in the interaction representation. This can be written 
in an alternative form, 

The proof of H25|l is obtained by setting F^, — F^o ^, so that 
dF, d , , , dF' , , ,dZi^ 

But by H19(l and the chain rule for the Poisson bracket, the final term can be written, 
{F^^, o Z;){Z^, GJ = {F', o Z„ GJ = {F,, GJ. (27) 
Rearranging the result gives 125|l . 

3. 5. Solutions for H^e and 

In view of (|25|l . the left hand sides of pt)a|l - pt)e|) can now be written DHne/De. 
In particular, ptiall is simply DHq,/ De — 0, which immediately gives _ffoe = C o Z,, 
where G is a function independent of e. Substituting e = 1 and the boundary condition 
shown in Tabled we find H ~ C o Z . Then substituting e = 0, we find G = Kq. In 
summary, 

iloe =KooZ,. (28) 

In particular, substituting e = 1 we obtain 

H = KooZ, Ko = Ho Z-^. (29) 

This completes the solution of and K to lowest order. 
The second order equation (|16&(l can now be written 

R^^l-{G,,Ho.h. (30) 

We use (|24f) in this, compose both sides with Z^^, integrate between e and 1, and use 
the boundary condition = at e = 1. The result is 

1 /"^ 

H2e = -^J de'{G,>,Hoe'}3oZ;;'oZ,. (31) 
Finally, setting e = 0, we have 

K2 - de{G,,Ho,}3 o Zr\ (32) 

Similarly, we solve the fourth order equation pCd) . finding 



K^ = -^l\e{G,,H2,hoZ-^ + j^l\e{G,,Ho,hoZ-\ 



(33) 



Clearly the solutions for Hnc and A'„ at any order n can be written in terms of integrals 
over lower order solutions. 

Let us now choose tl so that K will be a function of / at lowest order in h. We shall 
work backwards, first finding a canonical transformation Z such that Kq = H o Z~^ 
is a function of /. We then imbed this in a one parameter family , from which we 
compute Ge, Ge, Ue, and finally U. 
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3.6. Construction of Z via action-angle variables 

The desired canonical transformation Z can be specified in terms of the action-angle 
variables for the original Hamiltonian H and the harmonic oscillator. We let {A, (f) 
be the action- angle variables of H , according to the standard construction in classical 
mechanics, although we note that H may depend on h. The action is defined as a 
function of the energy by 

A{E) = 7^ / dpdx. (34) 
^ Ih<e 



271 



The integral is taken over the interior of the closed curve H — E {a. level set of H). 
The action vanishes at the fixed point, and is an increasing function of energy as we 
move away from it. Equation H34|l is the standard way to write the definition of the 
action, but, keeping in mind the warnings of [Appendix B} if we wish to think of A 
as a function on phase space, that is, a mapping : R'^ — > R then it makes no sense to 
write A{E). When we say that the action and energy are functions of one another, 
what we mean is 

H = fooA, (35) 

where fo : R ^ R is the function of a single variable expressing the relationship 
between energy and action. (The subscript will be explained below.) The function 
fo is invertible in the region of interest, so A = f^^ o H. Then 1)34(1 can be written 
more properly by picking a point z in the region in question, writing E = H{z), and 
then writing 

A{z) = {f^^oH){z)^f^\E), (36) 

instead of the left hand side of ((34|l . Having defined the action A, we then define 
the conjugate angle (j) by standard means in classical mechanics (through a generating 
function) . This involves a choice of origin (a point where (/> = 0) on each of the closed 
curves H = E in the region of phase space under consideration. This choice is smooth 
but arbitrary. We henceforth regard A and (p as specific functions : R^ R. 

Next we introduce the harmonic oscillator action-angle variables (1,9), where 
I is given by (jHJ and 9 is the conjugate angle (the geometrical polar angle in the 
phase plane, increasing in the clockwise direction). These are regarded as functions 
: R^ ^ R. Then the canonical transformation Z : R^ ^ R^ is defined by 

A = IoZ, 

37 

= 6* o Z, ^ ^ 

which uniquely defines Z, since {A, cj)), z — {x,p), and (/, 9) are invertible functions of 
each other (in the regions of interest). Then, using the definition (|29|l of Kq, we have 
Ko = HoZ-^ = fooAoZ-^ =foo I, (38) 

that is, Kq is the same function of the harmonic oscillator action / as iJ is of its own 
action A. This is reasonable, since canonical transformations preserve area, and the 
action is proportional to the area inside a level set. In particular, with this choice of 
Z, the level sets of Kq are circles. 

The smoothness of Z is relevant for the use of the Moyal star product series (|A.1|I 
which involves derivatives of functions on phase space. We believe that Z is smooth 
at all points of the relevant domain, under the assumption that H is smooth and 
has a generic fixed point. If the fixed point is not generic (for example, the quartic 
oscillator), then Z need not be smooth. These questions are discussed further in 
[Appendix C| 
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3. 7. Finding , and 

Now that we have Z, we imbed it in a smooth family Z^ with the boundary values 
shown in Table ^ Sjostrand and Zworski (2002) show that this can be done in a 
neighborhood of the fixed point, and Evans and Zworski (2004) give another proof 
that applies in the full domain. For later reference, we also define e-dependent versions 
of the action-angle variables, 

Ae^IoZ,, (39a) 

^eoZ„ (39&) 

with boundary values shown in Tabled Then we have 

iJoe - /O O A,. (40) 

All three Hamiltonians, iJ, iJoe and K are the same function (/o) of their own actions 
{A, Ae and /, respectively). 

Next we wish to find a function Ge such that (|19(l is satisfied for the given Z^. 
This can always be done, since that equation can be solved for the derivatives Gg^^, 
the components of a closed 1-form (hence exact, since the region is contractible). 
This is a standard result in classical mechanics (Arnold, 1989), which is summarized 
in component language in [Appendix E| The function is determined to within 
an e-dependent, additive constant. In the following we drop this constant, since its 
only effect is to introduce an e-dependent phase into Ue, which has no effect on the 
transformed Hamiltonian. 

Finally, given Gg, we transform it into the operator G^, and then define [7^ as the 
solution of 

subject to the initial condition [/g = 1 at e = 0. Then we set U = at e = 1. This 
completes the preparatory transformation (the construction of U such that K has a 
symbol that is a function of / at lowest order). We do not need to solve 141|l explicitly, 
since for the purposes of this paper we only need to calculate the effect on the symbol 
of a slowly varying operator when it is conjugated by U. But it is important to know 
that U exists, as we have shown. 

The preparatory transformation might have been carried out with oscillatory 
integrals coming from the integral representation of the Moyal star product, rather 
than in terms of a path through the group of canonical transformations. Indeed, 
we tried this approach initially, but found that it led to complicated algebra beyond 
lowest order that we were not able to organize to our satisfaction. Perhaps with more 
effort that approach could be cast into suitable form. 

The formalism we have presented is slightly simpler if we assume that the path 
through the group of canonical transformations, .Z^, < e < 1, is a one-parameter 
subgroup, that is, that G^ is independent of e. This, however, is a special assumption 
that we did not want to make. Moreover, the use of an arbitrary path allows us to 
study what happens when we vary the path, which leads to interesting conclusions 
(see below). 



4. Second stage transformations 



In the second stage (the second arrow in 0) we transform K into a new Hamiltonian 
M, such that the symbol M is formally a function of I to all orders in H. We do this 
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by Lie algebraic (power series) techniques that are similar to those used in classical 
perturbation theory (Dragt and Finn, 1976; Gary 1981), although here there are higher 
order Moyal brackets appearing as well as Poisson brackets. See also Littlejohn and 
Weigert (1993) for an example of a Moyal-based perturbation calculation applied to 
an adiabatic problem in quantum mechanics. 



4.I. The higher order transformations 

We apply a sequence of near-identity unitary transformations, each of which is 
responsible for making the symbol of the Hamiltonian a function of / at two successive 
orders of h. Only even powers of h occur in this process. The sequence is defined by 

m(") = k, 

M^*^ = M^^') Ul (42) 

etc., where 

;7„ = exp(-zn"-iG„), (43) 

and where Gn is the n-th order generator, assumed to have a symbol Gn that is slowly 
varying. Then, for example, the expression for M^^^ can be written as a series in h 
involving iterated commutators, 

M<^^^ =k-ih[G2,k]-^[G2,[G2,k]] + ..., (44) 

and similarly for Af'"*^ etc. Transcribing 144|l to symbols and using l)A.6(l . we have 

M(2) ^ K+h^G2,K}+h^ 2i^^'' + • • -(^^^ 

In a similar manner we write out commutator expansions for the higher order 
transformations in (|42|l . transcribe them into symbols, compose the transformations 
together, and substitute the expansion (|18() . We write the result in the form, 

M = Mo + h^M2 + h*lVU + . . . , (46) 

where M = M^°°\ the symbol of the final Hamiltonian after all the second stage 
unitary transformations have been carried out, and where 

Mo = Ko, (47a) 
M2 = K2 + {G2,Ko}, (476) 

M4 ^Ki + {G2, K2} - ^{G2, Ko}3 

+ ^{G2,{G2,Ko}} + {G4,Ko}, (47c) 

etc. Each M„ is slowly varying. 

We want M to be a function only of /. At lowest order we have this already. 

Mo = Ko ^ ho I ^ H o Z-\ (48) 

At second order, we wish to choose G2 in (|47bll so that M2 will be a function 
only of /, that is, independent of 9. In the next few steps it is convenient to bring 
back the abuse of notation rejected in [Appendix B| and to think of functions like K2, 
M2, etc. as functions of either z = {x,p) or of the action-angle coordinates (0,1), 
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as convenient. Then the Poisson bracket in ([47 6|) can be computed in action- angle 
variables, whereupon we have 

BC 

Ah^K^ + ^Loil), (49) 

where uj{I) — dKo/dl. Note that as a function, u; ~ Jq, since Kq — /o o /, so 
io{A) = dH/dA = /o(^)- Thus, U!{A) is the frequency of the classical oscillator with 
Hamiltonian H. If we now average both sides of (|49|l over the angle 9, we obtain 

A/2 = K2, (50) 

where the overbar represents the 9 average. The simple result is that M2 is just the 
average of K2, given by (|22Jl. 

Then subtracting H5U|) from (|49|) and rearranging, we obtain 

1 ^ ,,,, 

where the tilde represents the oscillatory part in of a function. Equation H51|) always 
has a solution G2 that is a periodic function of 9, that is, it is a single-valued function 
of (x,p), since K2 has a Fourier series in 9 without the constant term. Thus we have 
shown that it is possible to choose G2 in (j476|l such that M2 is independent of 9. 

The same structure persists at all higher orders. For example, taking the averaged 
and oscillatory parts of the fourth order equation l|47d) yields an expression for M4 
that is independent of 9 and a solvable equation for G4. This shows that it is possible 
to transform the original Hamiltonian H into a function of the harmonic oscillator / 
to all orders in h, at least in the sense of a formal power series for the symbol. 



4-2. Doing the e-integral 



The following steps require some notation and an important theorem regarding 
averaging operators that are explained in [Appendix D| The theorem in question is 
IjD.Sp . which we apply to using to obtain a useful form of the expression 

for M2: 

f 1 



M2 - / de{ (a 



1 

24 



de ( G 



oz: 



(52) 



where the diagrammatic notation is explained in [Appendix A[ The e-integration in 
([52|l can be done, yielding an expression independent of e, that is, independent of 
the path taken through the group of unitary or canonical transformations used in the 
preparatory transformation. 

First we transform the integrand of l|52|) as described in [Appendix F[ to obtain 



1 



M2 = — 



de 







d 



24 
1 d 
MdJ \luo I 



1 



uj o Af 



Hoe 



oZ: 



1 



de ( H, 



2 0e 



Hi 



Oe 



oz: 



(53) 
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where once we have transformed A^^ into / by composing with we can pull the 
factors depending on it out of the integral, since they are no longer e-dependent. Next 
we use the methods described in [Appendix G| to guess and prove that 

i-^ {Hoe ^Hoc) = {Hoe ^Ge^Hoe). (54) 
This makes the integral (|53|l easy to do, yielding. 




Let us call the two terms on the right hand side of H55|l the "H-teim" and the 
"-fTo-term." Since A'o = fo o I, the Moyal bracket in the /-iTo-term can be expanded 
out by the chain rule in terms of derivatives of /o and diagrams involving /. Wc find 

{Ko ^ Ko) = 2/^/^' (I^i^I) + (/ ^ /) , (56) 

where /q means /qo/, etc., and where some diagrams have vanished since (/ — /) = 
{/, /} = 0. The nonvanishing diagrams can be calculated using ((SJ, which gives 

{I^I^I) =21, (/=!:/) =2, (57) 

so the i^Q-term is a function only of / and the angle average in (|55|l does nothing 
to this term. Finally we take the /-derivative and compute the ii'o-tcrm explicitly, 
finding, 

SI/ £111 

Ko-tei-m^-f-if^I. (58) 

The intermediate Hamiltonian K is not unique, because of the choice of the path 
through the group of canonical transformations that connects the identity at e = 
and the given transformation Z at e = 1. More precisely, Kq = H o Z^^ is unique 
because it is expressed purely in terms of Z, but K2 and all higher order terms depend 
on Ze at intermediate values of e. Nevertheless, by (|50|l . if we vary the path Z^ while 
keeping the endpoints fixed, K2 can change by at most a function whose ^-average is 
zero, so that M2 remains invariant. Such a function can be written as the 0-derivative 
of some other function. These facts are proven in [Appendix H| 



5. The eigenvalues 

We have shown how to transform the original Hamiltonian H into a new Hamiltonian 
M whose symbol M is a function of / to any desired order in ft, and we have explicitly 
evaluated the first two terms Mq and M2 of the series for M. Let us write M„ = gn°I, 
thereby defining the functions (/„, so that M = g o I, where g = ga + ft^.92 + ^*.94 + ■ . •. 
In view of H48() we have go = fo, and 52 is given implicitly by H55() . 

As mentioned above, an operator is a function of / if and only if its symbol is a 
function of /. The two functions are the same at lowest order in h, but it turns out 
that they differ at higher order. These facts are proved in [Appendix I[ Thus, if we 
define a function / by M = /(/) and expand it according to / = + + fi'^fi + ■ ■ ■, 
then we will have fo = go but /2 ^ 52 • Thus fo defined this way is the same function 
introduced above in (|35|) . and we have M = fo{I) at lowest order. This is just what 
we guessed in the introduction, and it implies the usual Bohr-Sommerfeld formula, 
since the eigenvalues of / are (n + l/2)h. 
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5.1. The Bohr-Sommerfeld rule to higher order 

To carry the Bohr-Sommerfeld rule to higher order, it is necessary to find the relation 
between the symbol of an operator and the symbol of a function of that operator. 
This topic is discussed in [Appendix I| In the following we are interested in the case 
M = f{I) and M = g o /, so we will identify M and / with operators B and A of 
[Appendix I| respectively. Then (|I.4p gives the relation between functions / and g. 
Expanding / and g in even H series as above and using M„ = g„ o I, we can write (jl.4|l 
in the form, 

Mo + h^A'h + ... 



fir fill 
f2- — (I^I) (I^I-^I) 

16 ^ ^ 24 ^ ' 



+ 



= /o + ft^(/2-f -§)+..., (59) 

where /o means fo o /, etc., and where we use (|57|l . This implies Mq = /q o /, which 
we knew already, and allows us to solve for /2 by equating the final quantity in the 
parentheses with M2 in (|55|l . We see that the second order correction terms coming 
from (jl.4|) exactly cancel the i^o-term (|58|) . so that /2 o / is just the iJ-term of (|55|) . 



where we use (jP-Qjl . 

The eigenvalues of H are the same as the eigenvalues of M, which are given by 
/ o / evaluated at / = (n + l/2)h, or as we shall prefer to write it, f o A evaluated at 
A = {n + l/2)h (in this final step we are starting to confuse the functions /, A, with 
the values /, A). We compose foI = fQoI + h^f2 o I + . . . with Z and use 13 7|) . 
H48|l and H60|l to obtain (QJ, which is the Bohr-Sommerfeld formula including 0{h^) 
corrections. 

Equation is manifestly invariant under linear canonical transformations, since 
the matrix J'"' is invariant under conjugation by a symplectic matrix. Therefore, 
although this equation was derived in coordinates (a;,p) with balanced units of 
action"'^^^, the original units may be restored by a canonical scaling transformation, 
and the answer remains the same. 

In the case H = /2m + V{x), we have {H,H}2 = 2V"{x)/m, and agrees 
with the second order results of Bender and Orszag (1978), although we omit the 
details of the comparison. Equation ^ also agrees with the recent result of Colin de 
Verdiere (2004). We have also derived JQ) by a completely different method (a kind 
of WKB-Maslov method), and obtained the same answer. We beheve (Q) is correct. 

5.2. Action operators 

The formalism presented naturally suggests a definition of an "action operator." 
Let V be the overall unitary transformation resulting from the composition of the 
preparatory and second stage transformations, 

V = ...UiU2U, (61) 

so that 

M = VHV'^ = f{i). (62) 
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We then define an action operator 13 by 



(63) 



so that 



H 



V^f{i)V = f{B). 



(64) 



This is the relation whose expression in terms of symbols is the Bohr-Sommerfeld 
formula. It is straightforward to write out the symbol _B of i? in a power series in 
h. Our analysis of the multidimensional Bohr-Sommerfeld formula involves action 
operators in a more intimate way than the one-dimensional case. 

One can also transform creation and annihilation operators. Let a = (x -f 
ip)/{y/2h), = {x — ip)/{y/2h), so that / = {a^a + l/2)h, and define the unitarily 
equivalent operators b = V^aV, 6^ = V'^a'V. In this way many of the algebraic 
relations involving creation and annihilation operators for the harmonic oscillator go 
over to more general oscillators, for example, B = {b''b + l/2)h. 

6. Conclusions 

We conclude by presenting some comments on the present calculation. 

We could have expanded H in a power series in h, as in , and used the boundary 
conditions Hne — Hn at e = 1, which would have made all the symbols of this paper, 
Hne, Km M„, etc., independent of h. We did not do this because the odd powers of 
h in the expansion of H would complicate all subsequent formulas without otherwise 
raising any new, essential issues to be dealt with. The essence of the procedure we 
have given is one that operates only with even powers of h. 

In the calculation above there was a "miraculous" cancellation of the KQ-ievm. 
H58|l. where in one instance it arose as a consequence of doing the e-integral for K2, and 
in the second as a consequence of working out the symbol of a function of an operator. 
One suspects that this cannot be accidental. We will provide a deeper insight into 
this cancellation in our subsequent work on the multidimensional problem. 

The derivation of the multidimensional generalization of the Bohr-Sommerfeld 
formula (including order h? corrections), also known as the Einstcin-Brillouin-KcUer 
or torus quantization rule, requires new diagrammatic methods not considered in this 
paper. The answer is not an obvious generalization of the one-dimensional formula, 
and it involves some new geometrical issues for its interpretation. These topics will 
be the subject of a companion paper. 
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Appendix A. Notation for Moyal star product 

The Moyal star product A * B oi two symbols A, B is the symbol of the operator 
product AB. We write the h expansion of this product in the following notation, 



A*B 

^ — ' n\ 

n=0 

We call the bracket { , }„ that occurs in this series the "n-th order Moyal bracket" 
(other authors use this terminology to mean something else). This bracket is defined 
as follows. First, we define the Poisson tensor and its inverse by means of component 
matrices in the = {x,p) coordinates, 

Note that J^i, are the components of the symplectic form. We use J'"' or J^^ to raise 
and lower indices. This proceeds much as in metrical geometry, but one should note 
the sign change in 

X^Y^ = -X^Y^. (A.3) 

(In this paper we sum over repeated indices.) Next, we define 

{A,B}o=AB, (AAa) 

- A ^ J^"B (A.46) 

{A,B}2^A^^,J'^"rf'B,^p, (A.4c) 

{A, Bh = r'^r^r'^ (A.4d) 

etc., as required for IjA.ll) to be the standard Moyal series for the star product. Note 
that { , }i is the usual Poisson bracket. In this paper a bracket { , } without a 
subscript will be assumed to be a Poisson bracket. Note also that 

{A,B}„ = (-1)"{B,A}„. (A.5) 

Finally, note that if (7 = [A, B] , then the Moyal series for the symbol of the 
commutator is 

n=l,3,5,... ' ^ ^ 

= {{A, B} - ^{A, Bh + ^{A, S}5 - . . .) , (A.6) 

which defines the notation [A, S],. 

In this paper we make use of an alternative, diagrammatic notation for n-th order 
Moyal brackets and related expressions. For example, the ordinary Poisson bracket is 
be written 

{A,B}^ A^B, (A.7) 

where the arrow indicates differentiations applied to the operands A and B, connected 
by the J''" tensor. The base of the arrow is attached to the first index of J'"' and the 
tip to the second index. The operands can be placed in any position, as long as the 
arrow goes in the right direction: 

A B 

A-^B = \ = \ = B^A. (A.8) 
B A 
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But if the direction of the arrow is reversed, then there is a sign change, due to the 
antisymmetry of J^"^: 

A^B =-{A^B) , (A.9) 

which is the usual antisymmetry of the Poisson bracket. Similarly, the second Moyal 
bracket is given by 

{A,B}2= A^B = A^B . (A.IO) 

The two expressions on the right are equal because changing the direction of both 
arrows changes the sign twice. In this notation, the Jacobi identity is 

l(^A^B)^C]+[{B^C)^A]+[{C^A)^B] = 0, (A.ll) 

where the square brackets are only for clarity. The first term can be expanded out by 
the chain rule, which in diagrammatic notation gives 

(^A^B)^C = {A^B^C) + (C^A^B). (A.12) 

Similarly expanding the other two terms gives the vanishing sum of six diagrams, 
providing a diagrammatic proof of the Jacobi identity. 



Appendix B. Notation for functions 

In this paper it is convenient to use the (slightly nonstandard) notation / : A — > B to 
mean that the domain of function / is some suitably chosen subset of set A (in the 
standard notation, A itself is the domain). 

For the calculations of this paper it is important to avoid the usual abuse of 
notation in physics in which a function is confused with the value of a function. 
(Actually it is practically impossible to avoid this everywhere, but we shall do 
so wherever it is likely to cause confusion.) A "function" means a mapping, for 
example, H, H^, Ge, . . . : ^ M, and a canonical transformation is another mapping, 
Z„ Z, . . . : R2 ^ rpj^g components /x = 1, 2 of Z or will be denoted Zf" or Z^; 
each of these is a function : — » M. Functions will be denoted by bare symbols, 
H, Z^, etc., whereas values of functions will involve the specification of an argument, 
H{z), Z^*{zq), etc. It is also important to distinguish the identity map Id : ^ R^ 
from its value, which are the coordinates themselves. The identity map is defined by 

\d''{z) = z^. (B.l) 

One must also be careful about notation for derivatives. We use comma notation 

for derivatives since notation such as dA/dz^^ prejudices the choice of symbol to be 
used for the argument of the function. For example, the notation 

^Miz)) (B.2) 

is ambiguous; do we differentiate first and then substitute Z{z) for the argument, 
or substitute first and then differentiate? To avoid this problem, we write A,^ for 
the derivative of A, A,/^ oZ if we wish to differentiate first and then substitute, and 
(A o Z),^ if we wish to substitute first and then differentiate, where o represents the 
composition of two functions. The latter expression can be expanded by the chain 
rule, 

{AoZ),, = {A^oZ)Z';^. (B.3) 
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Poisson and Moyal brackets defined in (|A.4a|l - ljA.4d|) always denote functions. 
For example, the notation {A{z), B{z)} is meaningless, because it is only possible to 
take the Poisson bracket of functions, not numbers (the values of functions). On the 
other hand, {A,B}{z) is meaningful. 



Appendix C. The smoothness of Z 

The Moyal product rule IjA.ljl involves derivatives of symbols, and is not meaningful as 
it stands if the symbols are not smooth. We are assuming that H is smooth, but there 
is the question of the smoothness of the transformation Z as we have constructed it. 
We believe that under the our assumptions about H, the transformation Z is smooth, 
over a domain which is the open interior of a level set of H surrounding the fixed point 
but lying inside the first separatrix. We have not proved this, but in the following we 
present some considerations relevant to the question. We also present an example in 
which some of our assumptions about H are violated and Z is not smooth. 

There are two canonical transformations that are used in the construction of Z, 
one taking us from z = {x,p) to (0, A), and the other from z = {x,p) to {0,1)- 
The transformation Z defined by H37|l is the composition of one of these canonical 
transformations with the inverse of the other. These two canonical transformations 
are smooth except where 6 or <j) jumps from to 2-k, and except at the fixed point, 
where I = A — {) and or is undefined. Therefore Z is also smooth, except possibly 
at these places. 

Let the transformation from z = {x,p) to (0,1) be given by 

X — V2Ismd, 

(C.l) 

p — V 21 cos 9, 

which amounts to a convention for the origin of the angle 9 (it lies along the p-axis). 
This transformation is written without regard to the warnings of [Appendix B| but if 
we think of = {x,p) as values (g M) and 9 and / as functions, then the equation 
is put into proper notation by writing the left hand side as Id^(z) and 9 and / on 
the right hand side as 9{z) and I{z). Thus, lj('.l|l expresses the relation between the 
functions Id'^ and functions {9, 1). Now composing this with Z gives 



'2Acos( 

which gives an explicit representation of functions Z^ in terms of functions A and (f>. 
In particular, this shows that Z is continuous when jumps from to 2n, as long as 
A^O. 

The inverse transformation can be handled in a similar way. Let the 
transformation from z>^ = {x,p) to (0, A) be expanded in a Fourier series in 0, 

z^=^z,'^(A)e-^ (C.3) 

n 

where : M — > C are the expansion coefficients. This is subject to the same warnings 
about abuse of notation as (|C.ip . When these are straightened out and the result is 
composed with Z~^, we obtain 

(Z-^r = V(z^^o/)e»^ (C.4) 
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an explicit representation of Z^^. 

Now we wish to show that Z is smooth at the fixed point. This means that Z 
has derivatives of aU orders. In the following we present algorithms which we believe 
correctly give the derivatives of Z and of other functions at the fixed point, although 
we do not attempt to prove this in detail. 

Since H is smooth, it has an expansion about the fixed point, 

H{z) = ^Q^,z^z-' + ..., (C.5) 

where for simplicity we assume that the fixed point is at z = and that the constant 
term in the expansion vanishes, and where Q^i, is the positive definite Hessian matrix 
© evaluated at the fixed point. This series need not converge, but all coefficients in 
the series (the derivatives of H) are defined. It is convenient to manipulate such power 
series in a formal manner, since the rules for manipulating power series (multiplying, 
inverting, composing, etc.) are equivalent to the rules (chain, Leibnitz, etc.) for 
expressing the derivatives of new functions in terms of given derivatives of old ones. 

Birkhoff normal form theory (Birkhoff, 1927; Dragt and Finn, 1976; Eckhardt, 
1986) is a convenient way of developing a power series expansion of the transformation 
Z. This theory takes a Hamiltonian represented as a power series in z = (x,p), 
whose leading term is a harmonic oscillator, and transforms it into a function of 
the harmonic oscillator action /. The expansion (|C.5|I of H does not begin with a 
harmonic oscillator, but can be brought into this form by means of a linear canonical 
transformation L : R^, as shown by the theory of normal forms for quadratic 

Hamiltonians (Arnold, 1989). That is, there exists linear symplectic map L such that 

{HoL){z)^^ai{x^+p^) + ..., (C.6) 

where ai > 0, and where the ellipsis represents cubic and higher terms in a power 
series. In this step we rely on the positive definiteness of Q^i/. Then Birkhoff normal 
form theory provides another (nonlinear) canonical transformation : ^ R^, 
represented as a power series in z in which the leading (linear) term is the identity 
transformation, such that H o L o N is a function of /, 

H o Lo N = foo I = ail + a2l^ + . .. , (C.7) 

where /o is the same function introduced in (|34|l and H35(l . The coefficients ai, 02, 
etc. are determined by Birkhoff normal form theory, and they give the derivatives of 
function /o, which is smooth at the fixed point. We remark that in a one-dimensional 
problem such as this one, there are no resonance conditions so the Birkhoff algorithm 
can be carried to any order in the power series. 

It follows from if^Jl) and ((23 that A = loY, where Y = (Lo A)"\ so Y has the 
same effect on / as does Z in (|37|l . This does not mean that Y — Z, because Y does 
not necessarily satisfy the second of equations (|37|l . But 6* o y is an angle variable 
conjugate to A, so it differs from only by some phase shift 5 that depends on A. 
Assuming this phase shift is well behaved at A = (this is really an assumption of 
reasonableness on the definition of </)), there exists another canonical transformation 
S smooth at z = so that Z = Y o S satisfies both halves of H37|l . In this way all the 
derivatives of Z at z = may be computed. 

We note that if our conditions on H are not met, then Z need not be smooth at 
the fixed point. For example, the relation between action and energy for the quartic 
oscillator {V{x) = x*) is given hy H ^ cA^I'^ , so Ka = H o Z^^ = c'{x'^ +_p2)^/^ 
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where c and c' are constants. Thus, Kq is not smooth at the fixed point, and neither 
is Z. 

Appendix D. Notation for averaging operators 

This appendix develops abuse-free notation for the averaging operator introduced in 
subsection l4.1l Let Q : — > R be a function on phase space, treated as a Hamiltonian 
with evolution parameter a, 

^ = J^'-Q.^C^), (D.i) 

da 

where we assume Q is independent of a so the equations are autonomous (unlike 
the case of Ge considered above). Let : R^ — > be the associated flow, with 
components (Y^)''. The superscript Q indicates the Hamiltonian function generating 
the flow. The flow functions satisfy 

= jA"^(Q,,orQ) = {(rQf,g}. (d.2) 

da 

We will be interested in the case that Q is an action variable, /, A^, or A. 

For example, with Q = I, we have an advance map that advances the angle 
9 by a. That is, if a point z of phase space has action-angle coordinates {9, 1), then 
the point Y^(z) has coordinates {9 + a,I). Thus averaging over the angle 9 can be 
written as 

r^FoY^ = {F)o, (D.3) 



JO 2^ 

which defines the notation {F)g. Similarly, we define (F)^^ and (-F)^, using the 
advance maps Y^' and Y^. 

The advance maps and Ya are related by the following identity: 

Z^oY^' ^Y^oZ,. (D.4) 

In other words, angle evolution and e-evolution commute. We prove this by regarding 
both sides as functions of a at fixed e, and writing Xa and X'^ for the left and right 
hand sides, respectively. Note that Xa — X'^ at a = 0. The left hand side satisfies 
the differential equation, 

^ = (Z,. o r„^^)^^^ = (Z,. o Y^^){{Y^'-r. A.) 
da da 

= {Z,oY^%A,}^{Xa,A,}, (D.5) 

where we have used the chain rule property of the Poisson bracket. The right hand 
side satisfies 

dX' dY^ 

= -i^ o = {Fi, /} o = {X'a,A,}, (D.6) 
da da 

where we have used (|23|) and H39a|l . Since X^ and X'^ satisfy the same differential 
equation and the same initial conditions, they are equal, X^ — X'^, and the identity 
ljD.4|l is proven. It can also be written in the form, 

Y^' o Z-' = Z-' o ri. (D.7) 

It follows from this that for any function F on phase space, 

{FoZ-')^^{F)^^oZ-\ (D.8) 
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In view of Ij396|) this is a plausible identity. In particular, at e = 1 we have 

{FoZ-')^ = {F)^oZ-\ (D.9) 
To prove HU.8(I . we express the left hand side as an integral, and then apply (|D.7() : 



Appendix E. Function Ge exists 

Let Z^ be an e-dependent canonical transformation, defined on a contractible region. 
We wish to show that there exists a function such that H19|l is satisfied. Write 
S'^i, = Zj^j^ for the derivatives of Z^, which form a symplectic matrix. Then 

dZl 

or 



S^^r^G.^p, (E.l) 



dZ^^ dZ^ 

where we use the property of symplectic matrices, S^JS = J, where J is the matrix 
with components J^^. We must show that the second derivatives Gcp-y are symmetric. 
Differentiating, we find 

dZ^ dS^ 

Ge,i3-f = S" p^jJa^ +S°'/^Jafj, (E.3) 

The first term on the right hand side is symmetric in (/?, 7), since 

and the second term is also, as we see by differentiating S^JS = J with respect to e 
and juggling indices. Thus, the function Ge exists. 



Appendix F. A transformation of the integrand of ()52|) 

In this Appendix to save writing we drop the e subscripts on TJge, Ge, and 4>e, 
writing simply ifg, G, A and <j). The latter symbols, however, are not to be confused 
with the notation indicated in Table at e — 1 . We have placed this part of the 
calculation in an Appendix, to avoid confusion due to the notational change. 

Let us pick out the 0-average of the Moyal bracket in the integrand of H52|) and 
write it in an obvious notation, 

(c^iJo) - (G^Ho), (F.l) 

where the parentheses are only for clarity. We now introduce a technique for "breaking 
a bond" of an angle-averaged graph that is sometimes useful. The average of course 
depends only on A, if we think of it as a function of (0,^). We imagine evaluating 
this average at constant action A = a, which we enforce by inserting a (5-function and 
integrating over both A and (f>. This transforms (|F.1|I into 

m^SiA-a) (G^Ho), (F.2) 

where the integral is taken over a region of phase space that includes the level set 
A = a (an orbit of Hq). We then transform variables of integration to z = we 
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use dAd(f> = d^z (since the transformation is canonical), we write out one of the bonds 
explicitly, and we integrate by parts in the variable z'^: 



ZTT 



:i/o..) / 27 



+ 6{A-a) (G^f^.^Ho) 



(F.3) 



The second term in the final integral vanishes, due to the symmetry of G^ and the 
antisymmetry of J'^'^. In the first term we switch variables of integration back to 
(0, A) and do the A-integration, which gives 



27r dA 



G: 



Hn) =- — 



1 



Hn 



G: 



■ Hn 



where we use Hq = fo o A, that is, 14UII . and ui 
arrow in the final form. 



dA \u}o A 
/q, and change the direction of an 



,(F.4) 



Appendix G. An antiderivative for the integral ()53p 

In this Appendix we use the same notational simplifications as in [Appendix F| 

In guessing an antiderivative that will allow us to do the integral (|53|l . we must 
express the diagram Hq G =S Hq , which contains one G and three arrows, as 
D/De of some other diagram. We note by l|25|l that taking D/De of a diagram 
introduces both G and an extra arrow. Therefore taking the antiderivative must 
remove G and one arrow. The only diagram we can form from two copies of Hq and 
two arrows is Hq =^ Hq , so we compute, 

^^(i^o^i/o) ^(f^^i/o)-i(i/o^i^o)^G. (G.l) 
The first term can be written, 

(Ho^G)^Hq = {Hq^Hq- 

+ (Hq^G^ Hq) , (G.2) 

where we use (|16a|l and the chain rule, while in the second term of (|G.ip removing 
the parentheses provides a factor of 2, thereby cancelling the first term on the right 
hand side of (|G.2|) . As for the triangle diagram, it vanishes, as we note by writing, 





G 



Hq ^ — Hq 




(G.3) 



where in the first step we reflect about the vertical line and in the second reverse the 
directions of all three arrows. The overall result is (|54|l . 



Appendix H. The uniqueness of the intermediate Hamiltonian K 

In this appendix we study how the intermediate Hamiltonian K changes when the 
path Zf: through the space of canonical transformations is varied. To do this we 
compose with a near-identity, e-dependent canonical transformation that becomes 
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the identity at e — 0, 1. This is equivalent to replacing Zj^ with + SZ^, where 
6Z^ = {Z^,F^}, where is a small, e-dependent function such that = at 
e = 0, 1. The corresponding variation in the inverse function (Z~^)^ can be found by 
varying Z~^ o Z^ = Id, which gives 

5{Z-^r ^-J^'' F,^,oZ-\ (H.l) 

Then we vary H28|) to obtain, 

5Ho,^{Ho,,F,}. (H.2) 

Finally, to get (5Gc, we vary H19() to obtain, 

d 



de 



(5Zr) = {<5Zr,GJ + {Z^<5GJ. 



(H.3) 



The first term on the right hand side is {{Z^f^, F^},Ge}, while the left hand side is 



d 



" de 



Rearranging this and using the Jacobi identity gives 

0, 



where we use or, 



De 



De 



(H.4) 



(H.5) 



(H.6) 



where we drop a possible e-dependent constant. 

In the next few steps we adopt the same notational simplification mentioned at 
the beginning of [Appendix F| and in addition we drop the e subscript on i^^ and Z^. 
Then we combine P2ll . ()H.l|l ! (pL2)l and (pL6)l to obtain. 



SK2 = -- 



24 



de 



dF 



,Ho} -{{F,G},Ho}3 + {G,{Ho,F}h 

3 



de 

-{{G,Ho}3,F} 



o Z- 



(H.7) 



In this integral we perform an integration by parts, specified by 



D 
ITe 



{F,Ho}3 



dF 
Te 
d 
d'e 



,Ho} +{F,{Ho,G}h-{{F,Ho}3,G} 



oZ, 



{{F,Ho}3oZ-^) 



(H. 



which allows us to replace the first term of ljll.7|) with an exact e-derivative plus 
two more terms. The exact derivative can be integrated, giving zero because of the 
boundary conditions on F. What remains is 

r-l 



SK2 = - — 
24 



de 



{F, {Ho, G}}3 + {{F, Ho}3, G} ~ {{F, G}, Hojs 
+ {G, {Hq, F}}3 - {{G, Ho}3, F}] o Z-\ (H.9) 



We now use an identity related to the Jacobi identity for operators, itself a 
consequence of the associativity of operator multiplication. Let A, B and C be any 
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three operators, and write out the Jacobi identity [A, [B, C]] + cychc = in symbol 
form, expanding star commutators according to ljA.6|l . The leading order term is the 
Jacobi identity for the Poisson bracket, and the next correction term is 

{A, {B, Ch} + {A, {B, C}}3 + cyclic - 0. (H.IO) 

Using this in l)H.9|l allows us to write the integrand as 

{Ho, {G, Fh} o = [{lu o A){A, {G, F}^}] o 

= -(a;o/)A({G,F}3oZ-i), (H.ll) 

where we have expanded the Poisson bracket with Hq in action-angle variables. 
Finally, on restoring the e's we have 

SK2 ^ ^{u; o I)^ de{G,,F,hoZ-\ (H.12) 

The variation in K2 is an exact 0-derivative, as claimed, and M2 is invariant under 
variations in the path Z^. 

We do not know whether the space of symplectomorphisms we are considering is 
simply connected, but if not there arises the possibility of distinct paths that are 
not homotopic. Since M2 is unique, it must be that the difference in K2 along such 
paths is still an exact ^-derivative. 



Appendix I. Functions of operators vs. functions of symbols 



In this appendix we calculate the symbol of a function of an operator, in terms of the 
symbol of that operator, as a power series in H. We briefly describe a Green's function 
approach to this problem, which as far as we know was first presented by Voros (1977) 
and which is discussed further by Colin de Verdiere (2004) . In this appendix we adopt 
a general notation, in which A is any Hermitian operator, / is any function : R — *■ R, 
and B = f{A). The problem will be to find the symbol B in terms of the symbol A. 

Let a G C and let Ga = l/{a~A) be the Green's operator associated with A. The 
symbol Ga of Ga may be computed by demanding Ga * {a — A) = {a — A) * Ga — I, 
expanding Ga = Gao + h,Gai + fi^Ga2 + • ■ •, expanding the Moyal star product, and 
collecting things by orders in h. One finds that only even powers of h occur in the 
expansion of Ga, and that otherwise it is easy to solve for the leading terms. Through 
second order, the results are 



GaO 



Ga2 



1 



-A' 
1 1 



8a- A 
{A: 



1 



a- A' 
:A) 



.A 



A)3 



(A^A^A) 



(I. la) 



(1.16) 



This is a special case of the symbol of a function of an operator. For the general case, 
write B = f{A) in the form, 

da f{a) 



B 



Y^Tii a - A 



(1.2) 
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where the contour T runs from —oo to +00 just below the real axis, and then returns 
just above it. On taking symbols of both sides, this becomes 



The Green's function method becomes tedious at higher orders, but recently Gracia- 
Saz (2004) has found convenient methods for calculating the higher order terms, 
including the multidimensional case. It turns out that the fourth order term in ljl.4|l 
contains 13 diagrams. The Green's function derivation of (|1.4|l has required / to be 
analytic in a strip around the real axis, but Gracia-Saz has shown that the same 
expansion holds more generally. 

It was stated above that an operator is a function of / if and only if the symbol is a 
function of /. We prove this by noting that an operator is a function of / if and only if 
it commutes with the unitary operator U{t) — exp{~itl /H) for all t. This follows since 
the spectrum of / is nondegenerate. But the unitary operator U{t) is a metaplectic 
operator (Littlejohn, 1986), so when we conjugate an operator, A 1-^ U{t)AW (t), the 
symbol A is rotated in phase space. Therefore an operator commutes with all U if 
and only if its symbol is rotationally invariant in phase space, that is, is a function of 
I. 

The same thing can be proven at the level of h expansions. The general term of 
the series (11.411 involves diagrams composed of copies of A connected by arrows. But 
if A = I, then all diagrams with three or more arrows attached to any / vanish, since 
/ is a quadratic function of z. Therefore the only nonvanishing diagrams are linear 
ones and circular ones. A linear diagram with n Fs (two on the ends and n — 2 in 
the middle) vanishes if n is even, and is 2(— 1)^"~^^/^J if n is odd. A circular diagram 
with n Fs vanishes if n is odd, and is 2(-l)"/2 if n is even. Equation (|57|l is a special 
case of these rules. For now the point is that both these diagrams are functions of /. 
Thus the entire series (11.411 is a function of /, for any function /. 
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